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Abstract 

Starting from an asymptotically correct three-body Coulomb wave-function, we de- 
termine the effect of Coulomb final state interaction on the three-particle Bose- 
Einstein correlation function of similarly charged particles. We numerically estimate 
that the Riverside approximation is not precise enough to determine the three-body 
Coulomb correction factor in the correlation function, if the characteristic HBT ra- 
dius parameter is 5 - 10 fm, which is the range of interest in high-energy heavy ion 
physics. 
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1 Introduction 



One of the most important tasks of high energy heavy-ion studies is to prove 
the existence of the elusive quark-gluon plasma and to study the proper- 
ties of this predicted new state of matter [1]. Hanbury-Brown Twiss (HBT) 
interferometry[2-6] of identical particles has become an important tool as it 
can be used to measure the evolving geometry of the interaction region. The 
quantitative interpretation of the HBT-results depends however critically on 
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the understanding of the role of the Coulomb interaction between the detected 
particles as well as on the influence exerted on the selected system by the re- 
maining particles [7-9]. In the attempts to create the quark-gluon plasma we 
use higher and higher energies and larger and larger colliding nuclei. In these 
collisions the multiplicity of particles gets very large and the importance of the 
HBT-study of more than two particles will grow. Recently, much work [10-12] 
has been devoted to what can be learned from the system of three charged 
pions, however without paying sufficient attention on how the Coulomb inter- 
action might change the predictions. In the few existing experimental studies 
on the HBT effect in the three-charged pion system, see [13] and references 
given there, Coulomb effects have been taken into account in the so-called 
Riverside approximation [14], or sometimes in more elaborated versions of 
it [15]. However, these schematic conceptual ansaetze are acceptable only be- 
cause a more accurate multiparticle Coulomb interaction treatment has not 
been available. 

In this Letter we propose a proper treatment of the Coulomb interaction be- 
tween three charged particles. It is based on non-relativistic three-particle 
scattering theory, assuming that their velocities relative to each other are 
non-relativistic which, however, is just the kinematic situation where the HBT 
study of three particles is of interest. We, moreover, show that the Riverside 
approximation follows as a special case when the source of particle creation 
becomes very small, in the same way as the Gamow penetration factor is the 
limit for the Coulomb correction in the two-particle correlation function. An 
estimate of the magnitude of the corrections for source sizes to be expected 
at RHIC and LHC is presented using current data from NA44 experiment on 
S-Pb collisions at CERN SPS. 



2 Coulomb corrections to three-body correlations 

2.1 Three- charged particle wave function 

In order to treat correctly the Coulomb corrections to the three-body cor- 
relation function knowledge of the three-particle Coulomb wave function is 
required. As mentioned in the Introduction, we restrict ourselves to the case 
that the transverse momenta of the three particles in the final state are small 
enough to make a nonrelativistic approach applicable. Hence the problem con- 
sists in finding the solution of the three-charged particle Schrodinger equation 
when all three particles are in the continuum. 

Consider three distinguishable particles with masses m; and charges e^, i — 
1,2,3. Let Xj and k; denote the coordinate and momentum (three-)vectors, 
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respectively, of particle i. From these we construct in the usual manner the 
relative coordinate = x, — Xj and the relative momentum k.^ = (m^-kj — 
m,ikj)/(mi + rrij) between particles i and j, the corresponding reduced mass 
being /i^ = miirij/ (m, + rrij). 

The three-particle Schrodinger equation reads 



H + £ V i3 - E \ <i 2 k 3 (xi, x 2 , x 3 ) = 0, 



i<j=l 

where 



^ = E^>0 (2) 

j =1 

is the total kinetic energy for three particles in the continuum. H is the 
three-free particle Hamilton operator and 



Vy(ry) = ^(r y ) + ^(ry) (3) 

the interaction potential between particles i and j, consisting of a sum of a 
strong short-range (V^) plus the long-range Coulomb interaction (V^r^) = 

e i e j/ r ij)- 

The exact numerical solution of the Schrodinger equation (1) for E > is be- 
yond present means, partly for principal and partly for practical reasons. For 
a brief discussion of the related difficulties see [16]. But, at least in all asymp- 
totic regions of the three-particle configuration space is the form of its solution 
nowadays known analytically [17,18]. It is simplest in the asymptotic region 
usually denoted by fl and characterized by the fact that - roughly speaking - 
all three interparticle distances become uniformly large, i.e., r 12 ,r 2 3,r 31 — > oo 
(for a precise definition of the various asymptotic regions see [18]). 

A convenient and widely used representation of a three-charged particle wave 
function which coincides in Q with the correct asymptotic wave function is 
given by [19,20] 

<U( x i> x 2,*3) ~° < ( 2 +) (ri2)< ( 3 +) (r23)< ( 1 +) (r 3 i). (4) 

Here, (r^) is the continuum solution of the two-body Coulomb Schrodinger 
equation 
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-^ +i/ « (r « ) -5"l ,/ 'S +,(r «)= ' <5 » 



describing the relative motion of the two particles % and j with energy 
The explicit solution of (5) is known, 



^l +> (r tJ )=N i3 e l ^F[-t VlJ , 1; i(fcyr y - ky • r«)], (6) 

with iVy = e~ 7Tr,ij / 2 r(l+ir)ij), and 77^ = e '^ jAt ' J being the appropriate Coulomb 
(Sommerfeld) parameter. F[a, 6; x] is the confluent hypergeometric function 
and T(x) the Gamma function. When writing the asymptotic solution of (1) 
in the form (4), use has been made of the fact that in Qq, the short-range 
part of the two-particle interactions can be neglected and, hence, only the 
Coulomb potentials survive. 

Let us add a few comments. 

(i) It is to be emphasized that in using a wave function of the type (4), the 
three-body system is considered as a sum of three noninteracting two-body 
systems (each on its two-body energy shell). That is, neither correlations be- 
tween the motions of the three particle pairs nor off-the-two-body-shell effects 
are included, which obviously can be true at most for asymptotic particle sep- 
arations (in fact, it holds true in, and only in, Q ). 

(ii) The wave function (4) provides a well-defined prescription of how to go 
to arbitrary, in particular small, values of the relative coordinates. Though, 
it ceases to be solution of the Schrodinger equation (1) for non-asymptotic 
values of the relative coordinates; nevertheless, it has proved to be rather suc- 
cessful in describing (e, 2e) processes in atomic physics. But it should be kept 
in mind that the extrapolation out of the region flo implied by (4) is highly 
non-unique. In fact, several other ansaetze which, of course, coincide asymp- 
totically in fl with (4) have been, and are still being, developed. 

(iii) We mention that ansaetze for three-charged particle wave functions have 
been proposed which are correct in all asymptotic regions (see, e.g., [18]). As 
is to be expected, they contain correlations between the motions of the three 
particles. However, their form is rather more complicated and, thus, will not 
be used for the present investigation. For some more details and for further 
references see [21]. 

Hence, in the present investigation we assume the three-charged particle wave 
function to be given in all of three-body configuration space as 

^tU(xi,x 2 ,x 3 ) := ^< ( 2 +) (ri2)< ( 3 +) (r 2 3)< ( 1 +) (r3i)- (7) 
Here, Af is an (undetermined) overall normalization constant. 
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2.2 Specialization to three identical bosons 



Given a general three-particle wave function \& k * k2k (x]., x 2 , x 3 ) it is a simple 
task to specialize to various interesting situations. As a particular application 
we treat explicitly only the case of three identical bosons with unit charges of 
magnitude |e| and mass m. The properly symmetrized wave function is 

4t£k 3 ( X l> X 2, X 3 ) = ^= {<+ k2k3 ( Xl , X 2 , X 3 ) + <tk 2 k 3 ( X l, X 3, X 2 ) 

+^ktk 2 k 3 ( x 2> x i> x s) + ^ktk 2 k 3 ( X 2> X 3, xi) 

+ ^ktlk 3 ( X 3 ) X l ) X 2) + 4tk 2 k 3 ( X 3, X 2, X l)} • (8) 

By choosing for ^ k ^ k2k3 (xi, x 2 , x 3 ) the - at least asymptotically in Q correct 
- ansatz 



4tk 2 k 3 ( x i> x 2, x 3) = * k t k2k 3(x!,x 2 ,x 3 ), (9) 
we find for the approximate symmetrized wave function the form 



+ ^ 2 + WvCWvCW) 
+ < ( 2 +) (r 2 i)^ k 7 2 ( 3 +) (ri3)< ( ^ ) (r 32 ) 

+ ^ +) (r 2 3)< { 3 +) (r 3 i)^ ) (ri 2 ) 
+ ^ + V 3 i)< ( 3 +) (ri 2 )^ 3 ( 1 f) (r 23 ) 
+ <i +) (r 32 )< ( 3 +) (r 21 )^ ) (r 13 )}. 



(10) 



The case that all three particles are uncharged follows from (10) by substi- 
tuting for the Coulomb wave functions the plane waves e lkij ' riJ . We 
denote the corresponding symmetrized three-uncharged particle wave function 

by ^ESks^i.^.Xa): 



if/ ) 5 ( \ ^0 f Jk 12 -ri 2 ik 23 -r 23 ik 3 i-r 3 i , c iki 2 -ri 3 ik 23 -r 32 ik 3 i-r 21 

^kik 2 k 3 l x i> x 2> x 3j • — |e e e + e e e 

_|_giki 2 -r 2 igik 23 -ri 3 gik 3 i-r 32 _|_ g*ki 2 -r 23 g*k 23 -r 3 i g«k 3 i -ri 2 
_|_giki 2 -r 3 igik 23 -ri 2 gik 3 i-r2 3 _|_ giki 2 -r 32 ^ik 23 -r 2 i gik 3 i -ri 3 | (11) 
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where Ao is the appropriate normalisation constant. 



A simple approximation to (10) which incorporates at least part of the Coulomb 
effects is obtained by substituting for each of the two-particle Coulomb wave 
functions the term e lkij ' ri3 'N^, i.e., by neglecting the hypergeomet- 

ric function part in the exact solution (6) (this is justified if, for all pair of 
indices (ij) = (12), (23), (31), the arguments (A^r^- — • r^) of the hyperge- 
ometric functions are sufficiently small, as it happens, e.g., for small relative 
distances). This leads to 

^t£k 3 (xi,x 2 ,X3) » N U N 23 N 31 ¥^ 3 (^ 2 ^ 3 ). (12) 
For the probability density we arrive at 



C£k 3 (xi,x 2 ,x 3 ) _ « G 12 G 23 G 3l ^ 2k3 ( Xl ,x 2 ,x 3 ) " , (13) 

where the Gamov factors Gij = \Nij\ have been introduced. This is noth- 
ing but the well known Riverside approximation which consequently finds 
some theoretical justification within the framework of the exact nonrelativis- 
tic three-body scattering theory as a kind of lowest-order approximation. But 
this justification is rather weak. For, in order to "derive" the Riverside approx- 
imation, the wave function (7) or (10) is used in the region r 12 ,r 2 3,r 31 fa 
where it is certainly a much poorer approximation to the exact solution of (1) 
than at the larger values of the relative particle separations as applied in the 
following. Thus, a better justified expression is obtained by using the wave 
function (10). Nevertheless, also in that case an overall factor N 12 N 23 N 31 can 
be extracted from the r.h.s. of (10), which for the probability density yields 
again a factor G\ 2 G 23 G 3 i. Note that in the context of three-body scattering 
theory we do not find any support for the Coulomb correction factors proposed 
in [15]. 
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2.3 Application to high-energy heavy-ion collisions 



The correlation function measuring the enhanced probability for emission of 
three identical Bose-particles is given by: 



n n i i ^ iV 3 (ki,k 2 ,k 3 ) 

c ^' k3)= N^ms^m (14) 

where kj is the three-momentum of particle i. The three-particle momentum 
distribution is A^ 3 (ki, k 2 , k 3 ) and the single-particle momentum distribution 
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is denoted by JVj(kj). This correlation function is presently, due to meagre 
statistics, only measured as a function of the Lorentz invariant Q3, defined by 
the relation 

Ql = ^12 + ^23 + ^31 (^) 



where kij = ki — kj. In these equations is ki the four- momentum vector defined 
by the three-momentum kj and the mass of particle i, m,. 

We can now calculate the Coulomb effect on this three-particle correlation 
function using 



Kcouiomb{Qz) — z ~ J 2'(16) 



/ rf 3 Xip(xi)rf 3 X2p(x 2 )rf 3 X3p(x 3 ) 


^ktk!k 3 ( X l ; X 2, X 3) 2 


/ rf 3 X 1 p(x 1 )rf 3 X2p(x 2 )rf 3 X3p(x 3 ) 


*k?k 2 k 3 ( X l> X 2, x 3) 



where p{x.{) is the density distribution of the source for particle i, taken as 
a Gaussian distribution of width R in all three spatial directions, p(x) = 
(2^)3 exp[— (7^2 )]• This formulation makes it possible to estimate the Coulomb 
effect as a function of the radius parameter R on the three-particle correlation 
function, and to compare this estimate with that calculated by means of the 
Riverside approximation through K^ lomb (Q 3 ) = G12G23G31. 

To this purpose we use the NA44 data sample of three pion events produced 
in S-Pb collisions at CERN [22]. The pions of this sample are identified and 
well separated charged particles. All magnitudes of relative momenta |ky| are 
larger than 5 MeV due to experimental constraints and Q 3 is larger than 15 
MeV. 

We have calculated the Coulomb correction, i.e. KcouiombiQz)' f° r r& dius values 
R=5 and 10 fm, i.e. the range of interest for high energy heavy-ion physics, 
and compared the result to the Riverside approximation, see Figure 1. For a 
detailed discussion of these data see the NA44 publication [22] . 

We find that already for a source of radius 5 fm we need to take the de- 
tailed calculation into account as the difference to the Riverside approxima- 
tion is around 5-10%, furthermore the inset clarifies that this difference is 
(^-dependent, resulting in significant changes in any parameter extracted to 
characterize the correlation function. For a larger source radius, 10 fm, the 
difference between the estimates is even more pronounced. We have checked 
that in the limit R — > we recover the Riverside approximation. 

In Figure 2 we show the result of a corresponding calculation for the two- 
pion case [23] , illustrating that a finite size of the source reduces the Coulomb 
effect, as estimated by the Gamow penetration factor, much in the same way 
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as it reduces the Coulomb effect in the three-pion case from the Riverside 
approximation. 
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Fig. 1. The figures show the three-pion Coulomb correction factor KcouiombiQs) as 
well as the Riverside approximation. In the upper panel, the input radius was 5 fm, 
while in the lower, it was 10 fm. The insets display the ratio Coulomb wavefunction 
integration to the Riverside approximation. Lines are shown to guide the eye. 



3 Discussion: Possible areas of future application 



Finally we highlight some of the physical problems in high-energy physics 
where the above results may have a future application or where some mea- 
surements already indicate the possible influence of the many-body Coulomb 
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Fig. 2. The figures show the two-pion Coulomb correction factor -^cwomfc (^12) as 
well as the Gamow approximation. In the upper panel, the input radius was 5 fm, 
while in the lower, it was 10 fm. The insets display the ratio Coulomb wavefunction 
integration to the Gamow approximation. Lines are shown to guide the eye. 

distortions on the observables. 

1) In high-energy heavy ion physics nuclei are collided, that carry a large net 
positive charge and some of this charge may remain till the end of the reac- 
tion in the central collision zone. Such a residual positive charge is expected 
to distort the single-particle spectra of positively charged pions in different 
manner than that of the single-particle spectra of negatively charged pions. 
Namely, negatively charged pions are attracted to lower momenta, while pos- 
itively charged pions are repelled to higher values of the momenta. This effect 
was seen in the NA44 data in Ref.[24] and modeled in Ref. [9]. This problem 
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has a subtle three-body aspect, namely the correlations between the positively 
and the negatively charged pions induced by their Coulomb interaction which 
is also of long range, hence in principle cannot be neglected. In this case, we 
have a charge distribution (Z, +, — ) for which our wave function ansatz (7) 
can be utilized without the need of further symmetrization to determine more 
precisely the magnitude of the residual charge effects on the N[~\k) / N[ + \k) 
ratio of particle spectra. 

2) If a residual net positive charge plays a role on the distortion of charged 
particle spectra, its effects must also be seen on the like-charged two-particle 
correlation functions. For example, the two-particle intensity correlation func- 
tions are symmetric in the absence of a central residual charge, C < ^ ++ - ) = \ 
however, any net residual charge breaks the charge conjugation symmetry of 
the system, C^~ + ^ ^ \ This question could be properly studied also in 
the framework of a Coulomb three-body problem with charges (Z, +, +) and 
(Z, — , — ) with the two-body symmetrization of the Coulomb relative wave- 
functions of the two identical bosons in this system . 

3) The two-particle unlike-sign correlation function is also distorted by the 
Coulomb-effects of the residual charge. If we order the observed particles (pi- 
ons), e.g., according to the magnitude of their momenta, and fix this order 
during the determination of so that, e.g., the first particle is always the 
one with the smaller momentum, we predict that 

ci + '-\k u k 2 ) ^ct' +) O^iM) if Z res ^0 (17) 



This can be determined in a straigth-forward manner in experiments, either 
directly or by determining the charge-asymmetry of the two-particle distribu- 
tion functions, 

»(+-) = N^fakJ-Nt-^M) 

7V(+-)(k 1 ,k 2 )+7V(-+)( kl ,k 2 ) [6) 



Proper 3-body Coulomb calculations to estimate the magnitude of this effect 
can be performed with the help of the ansatz (7). 

4 ) Essentially, the method of 3) works also for the Coulomb correction due to a 
net residual charge for unlike-particle correlation functions. The study of these 
correlation functions was proposed recently by Lednicky and collaborators to 
access the temporal sequence of the particle emission [25]. This kind of physical 
information is not directly accessible to like-particle correlation measurements, 
and may provide a unique information on the dynamics of particle produc- 
tion. Hence it is of great importance to determine experimentally the value 
of the charge asymmetry parameter of the correlation functions, Eq. (18). If 
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this parameter is non-vanishing, three-body Coulomb-wave calculations have 
to be applied for unlike-particle correlations functions before determining pre- 
cisely the temporal sequence of particle emission from high-energy heavy ion 
experiments. 

5) In particle physics, Andersson and Ringner recently suggested to utilize 
the three-particle intensity correlation function C^ki, k2, k^) as a test of 
the hadronization mechanism to see the longitudinal stretching of the string 
field in the decay of quark-anti-quark jets[12]. The method to estimate the 
Coulomb effects imbedded in this problem is outlined in the present paper. 

Although this list can be continued we stop here because the above discussion 
is sufficient to show that the presented results are rather useful for a broad 
range of physical problems in high-energy particle and heavy ion physics. 



4 Summary 

On the basis of an explicit, analytically given form of the three-body Coulomb 
wave function that is correct in a large (asymptotic) region of three-body con- 
figuration space, we have developed a new method to systematically correct 
for explicit three-body Coulomb effects which is applicable to data analysis in 
a broad range of measurements in high-energy physics. The Riverside approx- 
imation has been established as a limiting case for vanishing source sizes. 

Specifically, we have worked out our approach for a system of three charged 
pions. There, we have numerically estimated that the Riverside approximation 
is not precise enough to determine the magnitude on the 5 -10 % level of 
the three-body Coulomb correction factor in the correlation function, if the 
characteristic HBT radius parameter is 5 - 10 fm, which is the range of interest 
in high-energy heavy ion physics. 
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